Enhancing Transport Efficiency by Hybrid Routing Strategy 
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Traffic is essential for many dynamic processes on real networks, such as internet and urban traffic 
systems. The transport efficiency of the traffic system can be improved by taking full advantage of 
the resources in the system. In this paper, we propose a dual-strategy routing model for network 
traffic system, to realize the plenary utility of the whole network. The packets are delivered according 
to different "efficient routing strategies" [Yan, et al, Phys. Rev. E 73, 046108 (2006)]. We introduce 
the accumulate rate of packets, rj to measure the performance of traffic system in the congested phase, 
and propose the so-called equivalent generation rate of packet to analyze the jamming processes. 
From analytical and numerical results, we find that, for suitable selection of strategies, the dual- 
strategy system performs better than the single-strategy system in a broad region of strategy mixing 
ratio. The analytical solution to the jamming processes is verified by estimating the number of 
jammed nodes, which coincides well with the result from simulation. 



I. INTRODUCTION 

Recently, the real transportation or communication 
systems such as the computer networks [l], H[ , power grid 
[3, 0|, airport line @, and so on, have attracted a lot 
of attention from scientists due to the discovery of the 
topological features of their self-induced structures. The 
complex network theory , as well as the tools inher- 
ited from nonequilibrium statistical physics [9| have been 
successfully applied to study the dynamical properties of 
these real systems. 

The common character for these transportation or 
communication systems is to perform certain functions 
by transferring objects among connected elements, which 
often take the form of large sparse network. Free traffic 
flow on these networks is key to their normal and efficient 
functioning. However, they may actually suffer from the 
overload or traffic jam, which always disable the system 
partially for a period of time, or even be fatal to the whole 
system due to the consequential onset of cascades of over- 
load failures Therefore, many recent studies on 
the traffic networks have analyzed the critical proper- 
ties of the jamming and congestion transitions [16h26| |. 
And, the schemes to promote the performance of traffic 
systems are chiefly from two aspects, designing efficient 
routing strategies (27rl3r| or, optimizing the topology of 
the underlying network [l9|, [37l - l40j . The objectives of 
these schemes are, on one hand, to avoid the onset of 
congestion and, on the other hand, to have short deliv- 
ery times. 

The routing algorithm proposed in recent works are re- 
lied on the structural properties, as well as the global or 
local information about the dynamical state of the com- 
munication networks [27l - f36| . For example, the works 
of biased random walk scheme introduce the probabil- 
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ity to visit node depending on its degree [28|, l29( , or the 
queue length of packets [3(| • The works of shortest-path 
scheme consider the paths with minimized distance from 
any pair of source and destination [3ll | . For this scheme, 
the central nodes (with highest connectivity) are highly 
overcongested, inducing the bottleneck of the communi- 
cation capacity. The expended version of the shortest- 
path scheme with "effective distance" involving the con- 
gestion state (queue length of routers) may bypass the 
congested nodes locally and thus improve the perfor- 
mance (32[. While, the work of efficient-path schemes 
[HI propose the routing table of paths with the mini- 
mum summary of fc' 3 , with a turnable /3. For the value 
of P = —1 this scheme can effectively redistribute the 
heavy load on central nodes to some of the lower-degree 
nodes, and the system can reach a more than ten times 
high capacity of that with shortest-path scheme. We can 
see that, for certain amount of traffic request, the way 
to promote the performance of the system is to take full 
advantage of all kinds of nodes. 

These aforesaid researches, have discussed the system 
with pure routing strategy. While, how the diversity of 
routing strategy performs is really of curious, and the en- 
hancement of transport capacity by better exertion of all 
nodes in the system might be expected. In this paper, 
we put forward a mechanism that the communication 
system possesses of two different routing strategies. Here 
we make use of the simple fixed routing scheme, i.e., the 
efficient-path schemes proposed in Ref . [33| , and consider 
the routing strategies to be denoted by different /3. Then, 
the transport system with this multi-strategy protocol 
will send packets according to different fixed routing ta- 
bles of efficient-path schemes. Though the fixed routing 
algorithm becomes impractical in huge communication 
systems, it is still widely used in medium-sized or small 
systems [4l|, |42[ , for its obvious advantages in economical 
and technical costs, compared with the dynamical rout- 
ing algorithm and information feedback mechanism. In 
this case, the diversity of the fixed routing strategy is, of 
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course, practical if it performs better than pure-strategy 
system. Actually, through our study, we see that the 
multi-strategy system may perform better than that of 
the pure strategy system. 



II. TRAFFIC MODEL 

In our traffic model of dual-strategy routing protocol, 
the packets with given sources and destinations will be 
sent according to two different fixed routing tables of 
efficient-path schemes (EPS). For the EPS proposed in 
Ref . [33l] . node i in the graph are weighted by to* = fcf . 
hi is the degree of node i, and /3 can be considered as 
the label of "routing strategy" . A packet with source j\ 
and destination j 2 will choose a minimum sum of weight, 



SiGo- , route in the graph. (Tj 1 j 2 is the path from j\ 
to ]2- Adjusted by the parameter /3, the single-strategy 
system will partial to certain kind of nodes in routing, 
and may also leave some space to improve the perfor- 
mance further. In our dual-strategy model with two 
strategies /3i, and 02, packets are assigned to the two 
corresponding routing tables, with probability 1 — p and 
p, respectively. Here we name p as the mixing rate. Here, 
for p = (or 1), the system returns to the single-strategy 
system with /3 = fl\ (or /3 2 ). 

Similar to the former work, at each time step, R pack- 
ets enter the system with randomly chosen sources and 
destinations. The delivery capacity of each node is C, 
and we set C — 1.0 for simplicity. The maximal queue 
length of each node is assumed to be unlimited, and the 
first-in-fist-out discipline is applied at each queue. Once 
a packet come to its destination, it is removed from the 
system. 

In the previous study, the phase transition of traffic 
flow is described by the the order parameter [iff , 
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where AW = W (t + At) — W (t) , with (■) indicating aver- 
age over time windows of width At, and W(t) is the total 
number of packets in the network at time t. The criti- 
cal value R c (the packet generation rate) where a phase 
transition takes place from free flow to congested traffic, 
can reflect the maximum capability of a system. 

The behavior of the critical point R c on different net- 
works can be simply explained by their different between- 
ness centralities (BC) distributions dHHEl]. The BC 
of a node i for the single-strategy EPS system [33| is 
defined as, 
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where &j 1 j 2 {/3) is the number of routes going from j\ to 
j2, according to the EPS routing table with /3; While, 
(Jj 1 j 2 (P, i) is the number of those also passing through i. 




FIG. 1: The accumulate rate n as a function of /3 in the 
single-strategy system, for the systems of R — 50 and 60. 
The results are averaged over 10 realizations for 20 networks, 
with size N = 1225. 



The critical value R c can be estimated by the maximal 
BC as, 



Rc 



C ■ N ■ (N - 1) 

Max[ gi ((3)] 



(3) 



where Max[gi(f3)] is the maximal BC of the system with 
strategy (3. 

For the dual-strategy system with strategies j3\, P2, 
and probability p, the efficient BC of one given node i is, 

G0 u p 2 ,p) = {l-p)-gM+p-9M (4) 

Then, we have the load of node i, assigned from the whole 
transport requirement of the system as, 



U = 



G t (fii,p2,p)-R 
N ■ (N - 1) 



(5) 



The load of node increases as the R is increased. There- 
fore, the critical value R c can be estimated as, 



Rr.= 



C-N-(N-l) 

Max[G,(/3i,/32,p)] : 



(0) 



here, Max[Gi(j3i, 02, p)) is the maximal efficient BC of 
the dual-strategy system. 



III. SIMULATION RESULT AND ANALYSIS 

The communication networks typically show a scale- 
free (SF) distribution for the number of links departing 
from and arriving to a system element. In this paper, 
we choose Barabasi-Albert (BA) network as the commu- 
nicating network (45j . For this network model, starting 
from mo = 3 fully connected nodes, new node with m — 2 
is added in the existing network in turn, until the network 
size N = 1225. The network average degree (k) = 4. 
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For the single-strategy system, the phase transition 
from free flow to congested traffic has been discussed [33[ . 
When the value of R increases over R c , the number of ac- 
cumulated packets get to increase with time (i.e., a phase 
transition takes place from free flow to congested traffic). 
Similarly, for the multi-strategy system, the phase tran- 
sition also takes place. The effect of different strategies, 
in free flow phase, is merely inducing the difference of 
packet deliver time. While in the congested phase, much 
more diversified phenomenon appears. We mainly focus 
on the congested phase as follows. 

Firstly, let us revisit the behavior of the single-strategy 
system in the congested phase. According to the work of 
Yan (33|, the largest R c (around 43), i.e. the best perfor- 
mance of the system, is achieved with strategy ft = 1.0 on 
BA network of N — 1225 and (k) — 4. From systematic 
simulation of various (3 systems in congested phase, we 
notice that the number of accumulated packets increases 
linearly with t. Namely, the accumulate rate rj is a con- 
stant (with small fluctuation). In Fig. Q] we shows rj as a 
function of strategy /3, with R in the region of congested 
phase (R = 50 and 60, larger than R c ). It is necessary to 
emphasis that, although the so-called congestion occurs, 
there still are, on average, R — rj packets successfully de- 
livered to their destinations per unit time. This number 
is actually much larger than R — R c . That is to say, while 
some nodes are jammed as R > R Cl a noticeable part of 
transport function still holds in the system. This actu- 
ally is realized from two aspects, (1) the "free flow" still 
takes place on the paths which are not entangled with the 
jammed nodes, and, (2) the packets through the jammed 
nodes are not stopped but just delayed. 

We may say that, the parameter R c merely distin- 
guishes the so-called free and congested phases, which 
actually indicates the free or jammed state of the most 
"fragile" node [see Eq. ((3|)]. R c can not reflect the ex- 
tent of congestion, and the impact of the jammed nodes 
to the performance of the system. However, the accumu- 
late rate, defined as, 



?y = lim 



AW 



t^oo At 



(7) 



is a good parameter to measure the performance of the 
system in the congested phase. The smaller r\ denotes 
better performance of the system, rj is the sum of in- 
dividuals' 77- over the whole system as, -q — YliVi^iVi)' 
Here, H(-) is the Heaviside function, and 77- is the indi- 
vidual accumulate rate of node i, namely, the increase 
rate of the queue length of packets at node i per time 
step. 

From Eq. ©, we can get the analytically expression 
of r)l as, 
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C = 



Gj((31f32,p) ■ R 
N ■ (N — 1) 
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(8) 



with Li the load of node i assigned from the whole trans- 
port requirement. We may notice that as R is increased, 
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FIG. 2: (Color online.) The accumulate rate rj for the dual- 
strategy system as a function of mixing ratio p of the two 
strategies /3i and fa- Here, in (a), P2 is fixed to be 1.5, and in 
(b), /3i is fixed to be 0.5. The rj of the single-strategy system 
with optimal ft = 0.9 (the red dot line) is also plotted for com- 
parison. The results shown are averaged over 10 realizations 
for 20 networks, with size N — 1225, and R = 60. 



Li may increases over the capability C and thus r}[ in- 
creases from negative to positive. 

In Fig. [TJ the non-monotonic behavior of rj implies 
that the medium / 3 sy stem performs better, similar to 
the results in Ref. [33| from the relationship between R c 
and /3. 

Then, we will analysis the behavior of the dual-strategy 
system with /?i and /?2 in the congested phase. The pack- 
ets are assigned to the two strategies with probability 
1 — p and p, respectively. Figure [5] plots rj of the system 
as a function of p. Here, for p — (or 1), rj returns to that 
of the single-strategy system with ft = fli (or /3 2 ). We can 
see that, the mix of different strategies is nontrivial and 
of interest. Take the system with j3i = 0.5 and /?2 = 1.5 
in Fig. [DJa) as an example, for certain medium value 
of p, it performs even better than the optimal state the 
single-strategy system achieves with f3 = 0.9 (which is 
also plotted by the red dot line in Fig. [5]). Furthermore, 
as has been shown in Fig. [21 it is also noteworthy that, 
when fix and fit are chosen from each side of 0.9, there 
always exists an optimal configuration p, which performs 
better both than the single-strategy systems of /3i and 
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This can be understand as follows. To design routing 
strategy for the network transportation, there are two 
factors that should be considered. (1) To bypass the hub 
nodes which are obviously of heavy burden and prone to 
jamming. (2) To choose shorter path to reducing deliver 
time, which is conducive to reduce the occupation (life 
time) of packets to the resources and thus avoid jam. The 
system deliver efficiency can be improved from the trade- 
off of these two factors. However, they are inconsistent 
in the communicating network with heterogeneous topol- 
ogy. Take the single-strategy system in congested phase 
as an example (see Fig. Q] the curve with R = 60), as 
fi is increased from 0, the traffic through the hub nodes 
are bypassed to the other smaller degree nodes, while the 
lcngthes of the pathes adopted are prolonged, which in- 
creases the probability of jamming for the other nodes. 
The system with /3o around 0.9, to certain extent, is com- 
patible of these two factors, and thus achieves the optimal 
performance. As /3 is increased further, the utility of the 
hubs is not sufficient, while the left parts of the system 
are overworked. Actually, To take a full advantage of 
each node in the system will return better performance. 
Therefore, for the dual-strategy system, the strategy in- 
clined to the hubs ((3 < (3q) and that inclined to the 
small nodes (/3 > /?o) may complement each other and 
perform better than the single strategy one. Thus non- 
monotonous rj can be observed when the (3 from both 
side of /3o are mixed. 

The effect of multiple strategies in the congested phase 
can also be understand analytically from the so-called 
equivalent generation rate. In this routing strategy, pack- 
ets at the head of the queue on node i will be delivered to 
the next node j according to the routing table, no mat- 
ter node j is idle or jammed. Current server also has this 
properties. In this case, congestion in the system will 
not spread out. Furthermore, counterintuitive, conges- 
tion will make the system more "empty" . In each time 
step, rj more packets will queue at the jammed nodes, 
and as a consequence, the load of the other nodes will 
be lighten, as if the generation rate for the subsystem of 
these nodes is reduced to a smaller one R* , which we 
name as the equivalent generation rate. Here, we have 



R* 



R-rj. 



(9) 



Different from the case that the servers abandon packets 
when the queue length is over a threshold, in our model, 
the queuing packets are not abandoned, and will finally 
be send to their destination. 

We sort nodes by the values of their individual accu- 
mulate rates in descending order, as r/^ > i] 2 > ■■■ > r/' N . 
FromEq. ©, we know that, when R is increased from 0, 
all these 77^ increases from — C. As soon as the maximum 
one, r/'i, increases from negative to positive, the system 
transform from free phase to congested phase. Suppose 
that 772 < 0, there are n[ packets detained at the 1st node 
per time step. Then, the equivalent generation rate for 
the subsystem (exclude the 1st node) is R* = R — r][. 
As R is increased further, the left nodes will be jammed 
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FIG. 3: The number of jammed nodes from analytical and 
simulation results, for the dual-strategy system with /3i = 
1.5, @2 = 0.5. The sample data of analytical results (red open 
circle) are from 10 different networks, and that of simulation 
results (black open square) are from 50 realizations of traffic 
on these 10 networks. The average number of analytical and 
simulation results (red solid circle and black solid square) are 
averaged over the corresponding sample data. The system is 
of size N = 1225, and R = 60. The analytical and simula- 
tion results (the red and black dot lines) from single-strategy 
system with ft = 0.9 are also plotted for comparison. 



one after another (i.e., have positive T}[). Accordingly, we 
may propose the theory to predict the number of jammed 
nodes, and the accumulate rate of the system 77 from two 
perspective. 

On one hand, from Eqs. (J5J) and ([9]), we get, 



R* = R-J2l 



- L N-(N-l) 
with the following constraint applies: 

Gi(/3i,/3 2 ,p) ■ R* ^ n 
Ll N-(N-1) > °' 

G I+ ^ u p 2 ,p)-R* 
A - (TV - 1) 



C] (10) 



fill 



L1+1 — 



< C. (12) 



By solving this problem, we can get the number of 
jammed nodes /, and 77, for given values of R, /3i, f3 2 
and p. 

On the other hand, we focus on the detailed process of 
successional jamming which gradually modifies the equiv- 
alent generation rate R* , as well as the load L.; of the left 
nodes. The iterative procedure of R* can be written as, 



Rl = R- 



Gi(/3i,/3 2 ,p)-fl 
A - (TV - 1) 



C 



(13) 
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The iterative formula is, 

Ri = Ri-i n-(N-i) — ' ' = ' ' ' '"' 

(14) 

i?* and L\ decrease as the nodes of large load is jammed 
one after another, until 



L', 



G 7 (/3i,/3 2 , P )-^_ 
N-(N-l) 
G/+i(/3i,/3 2 ,p) 



>C, 



N ■ (N - 1) 



< C, 



(15) 
(16) 



Different from Eqs. (HHJ) to (H2]), Eqs. flU]) to (H]) depicts 
that the jamming of the first I nodes steps down R* 
gradually until the value R}, where the (/ + l) th node, 
as well as all its following nodes, is capable of treating 
with its load. Here, from the perspective of successional 
jamming process described by Eq. (fl4|) . one can also get 
the number of jammed nodes /, and t], analytically. 

In Fig. [31 we plot the analytical and simulation results 
of the number of jammed nodes / in the dual-strategy 
system with /3i = 1.5 and — 0.5. It can be seen 
that, the average number of jammed nodes from analy- 
sis (red solid circle) coincides well with that from sim- 
ulation (black solid square). Interestingly, the value of 
/ also behaves non-monotonically and achieve the mini- 
mum around p — 0.5, which is similar to the accumulate 
rate rj of the same system shown in Fig. [2] Additionally, 
the analytical results from Eq. (TlOl) and Eq. f| 14[) are 
very close to each other, thus in Fig. [3] we merely plot 
the results from Eq. (TT41 . 

Here, we can also understand the non-monotonic be- 
havior of / from the following perspective. The packet 
generation rate R can be divided into two parts, the pack- 
ets using routing table of (3i is R& 1 = (1 — p)R, and that 
of /?2 is R^ 2 = pR. From Eq. ([5]), we can get the corre- 
sponding loads of node i from these two parts of packets, 
denoted by if 1 and if 2 (with U = if 1 + L? 2 ). For the 
case that the mixing rate p = 0, we have R& 1 = R, and 
the jamming of nodes are all ascribed to the queue of /3i 
packets. As p is increased from 0, the RP 1 , as well as 
the if 1 decreases, while that of /?2 increases. If the /?2 
packets prefer to use those complementary nodes instead 



of the nodes already jammed by 02 packets, the number 
of jammed nodes / will decreases with p. However, as p 
is large enough, the increase of load Lf 2 from /?2 packets 
induces new jamming of nodes. Therefore, we can see the 
non-monotonic behavior of the number of jammed node, 
when the dual-strategy system is composed of the two 
strategies from either side of (3q. 



IV. CONCLUSION 

In summary, we propose a hybrid routing strategy for 
the networked traffic system, which is proved to be a 
doable and effective way to enhance transport efficiency. 
Compared with the efficient routing strategy [33j, the 
hybrid routing strategy can make better use of the re- 
sources in the traffic system, while there appears no in- 
crease in its algorithmic complexity. The performance 
of the dual-strategy system can be optimized by mod- 
ulating the mixing rate of the packets, in case that the 
two strategies share fewer key nodes. Here, we introduce 
the accumulate rate rj to denote the performance of the 
communication system in congestion phase, which shows 
richer phenomena than the critical generation rate R c . 
Furthermore, we get analytical descriptions to the jam- 
ming processes by the accumulate rate r\ and the equiva- 
lent generation rate R* . The number of jammed nodes es- 
timated from analytical formula coincides well with that 
from simulation. 

While our model is based on computer networks, we 
expect it to be relevant to other practical transport 
processes in general. Actually, in real system, the hybrid 
routing is worthy of considering, for the reason that 
the sources and characters of massages delivering or 
spreading in complex systems are diversified, which 
induces the hybrid of various transportation modes. In 
view of the common features for the networked traffic 
and spreading, our work may shed some light on the 
research of packet delivery in technical networks, as well 
as the rumor and opinion dynamics in social networks. 

We gratefully acknowledge T. Zhou and X. Li for help- 
ful discussions. 
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